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ABSTRACT 

We examine properties of thermal radiation emitted by strongly magnetized 
neutron stars (NSs). In particular, we show that the pulsation amplitudes of the 
energy-integrated flux are an increasing function of the absorption column density to 
the source. This is especially important for the interpretation of the Anomalous X-ray 
Pulsars (AXPs) as cooling neutron stars with high magnetic fields. We show that the 
high-pulsation amplitudes observed in these objects are consistent with cooling models, 
if the large amount of absorption to these sources is taken into account. We also show 
that cooling models imply inferred radii of the emitting regions on the order of ~ 5 — 6 
times smaller than the actual NS radii, again in agreement with observations. 

Subject headings: stars: neutron — X-rays: stars 

1. Introduction 

In the past few years, X-ray observations have revealed a new class of X-ray sources, the 
so-called anomalous X-ray pulsars (AXPs) (Mereghetti &: Stella 1995; van Paradijs, Taam & 
van den Heuvel 1995). AXPs are sources of pulsed X-ray emission, with persistent luminosities 
L x ~ 10 35 — 10 36 erg s _1 and soft spectra. Their periods lie in a very narrow range, between 6 and 
12 seconds, and their characteristic ages are of order 10 3 — 10 5 years. Their distances are usually 
very large, as inferred from the typically large column densities to them (Ah ~ 10 22 cm -2 ). 

Due to the long period of these objects, it is clear that their X-ray luminosity cannot be 
rotation-powered. Two broad classes of models have emerged to explain their X-ray emission. In 
one, the radiation is assumed to be powered by accretion from binary companions of very low mass 
(Mereghetti & Stella 1995), from the debris of a disrupted high mass companion (van Paradijs, 
Taam & van den Heuvel 1995; Ghosh, Angelini & White 1997), from the interstellar medium 
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(Wang 1997), or from the debris falling back after the supernova explosion (Chatterjee, Hernquist 
& Narayan 2000; Alpar 1999, 2000; Marsden et al. 1999). In the second class of models, AXPs 
are hypothesized to be isolated, ultramagnetized neutron stars, with field strengths in the range 
10 14 — 10 15 G, a.k.a. "magnetars" (Duncan & Thompson 1992). The X-ray luminosity could then 
be powered by either magnetic field decay (Thompson & Duncan 1996; Heyl & Kulkarni 1998) or 
by residual thermal energy (Heyl & Hernquist 1997a,b). 

There are several observational constraints that can be used to discriminate among various 
models. The presence of an optical counterpart and observable Doppler shifts in the frequency 
of the X-ray pulses would suggest the presence of an optical counterpart, but this has not been 
observed so far. Optical and longer wavelength emission (Hulleman et al 2000; Perna, Hernquist & 
Narayan 2000) would indicate the presence of an accreting fossil disk. Any viable model, however, 
must also be consistent with the observed high pulse fractions and small emitting areas inferred 
for the X-ray radiation from these objects. 

In this Letter we examine whether the observed pulse fractions and inferred emitting areas are 
consistent with a model for AXPs in which their X-ray emission is powered by thermal radiation 
from highly magnetized neutron stars. We show that the amplitudes of the pulsations for the total 
flux observed in a given band are strongly sensitive to the amount of absorption to the source, 
and, more specifically, increase with it. To understand this point, one first needs to notice that, 
for typical local blackbody (BB) emission, the energy-dependent pulse fractions increase with 
energy, as a consequence of the fact that the hardness of the blackbody spectrum becomes larger 
at higher energies. This point was well elucidated by Page (1995) and Page & Sarmiento (1996): 
the ratio of the BB fluxes emitted at energy E by two regions of areas A\ and A 2 with respective 
temperatures T\ and T 2 is 

F BB (£, Ji) _ A 1 exp(£ , /fcT 2 ) - 1 

F BB (E,T 2 ) A 2 exp(£7/fcTi) - 1 ' U 
which is an increasing function of E if T\ > T 2 . This property remains valid for atmospheric 
models with beaming oc (cos#) n , and, furthermore, the increase of the pulse fractions with energy 
becomes stronger with increasing n. 

Absorption by the intervening interstellar medium (ISM) does not affect the energy dependent 
pulse fractions (though this is rigorously true only for a perfect detector), as photons emitted 
at a given energy are reduced by the same amount independently of the phase of the rotation. 
However, pulse fractions measured from flux which is integrated over a given band do depend 
on absorption. Absorption, in fact, reduces the low-energy part of the flux with respect to the 
high-energy one. As a result, the main contribution to the pulse fractions is weighed more and 
more by the high energy component as the amount of absorption is increased. 

We compute the expected pulse fractions and inferred areas for typical parameters of AXPs 
and a thermal cooling model for their emission, which includes processing through an atmosphere. 
We show that the inferred properties are consistent with observations. 
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2. Model 

Our model is that of a neutron star (NS) cooling through an accreted envelope and with an 
intense, dipolar magnetic field. Heyl & Hernquist (1998b) showed that, for B p ^ 10 14 G, the flux 
transmitted through the envelope is well described by F oc cos 2 ip, where ip is the angle between 
the radial direction and the magnetic field. For a dipolar field, cos 2 ';/' = 4 cos 2 # p /(3cos 2 9 P + 1) 
(Greenstein & Hartke 1983), with 9 p being the magnetic colatitude (angle between radial direction 
and magnetic pole). It is given by 

cos 6 P = cos 6 cos a + sin 6 sin a cos <fi , (2) 

where a is the angle that the magnetic pole makes with the line of sight, and (6, (ft) are the 
usual colatitude and longitude in spherical coordinates. In our calculation, we assume a = 90° 
(orthogonal rotator). Therefore our results on the pulsation amplitudes must be interpretated as 
upper limits. We assume that the local emission, n(E,T), is a blackbody spectrum modified by 
the presence of an atmosphere, for which we adopt the semianalytical model of Heyl & Hernquist 
(1998a), but with the inclusion of limb darkening 0. This leads to a beaming of the radiation 
oc cos 5, with S being the angle between the normal to the NS surface and the direction of the 
photon trajectory. We will also consider cases with a more intense beaming oc (cos<5) n , but not 
extreme (i.e. we restrict our attention to n < 2). Static general relativistic effects (i.e. light 
deflection and gravitational redshift) are fully taken into account in our calculations. 

Let us define e _As = y/1 — R/R s , where R is the radius of the NS star, and R s = 2GM/c 2 
its Schwarzschild radius. Here we assume M = 1.4Mq, and consider the range of radii 
2 < (R/R s ) < 4, all compatible with the currently available models for the NS equation of state. 
Let D be the distance from the star to the observer, and iVn the intervening column density. The 
flux measured by an observer at infinity is then given by 

m m = =%g^e-<^» /; 2^ f I ,„(., « n\Ee^ ; T.(», «] , (3) 

in units of phot cm" 2 s" 1 keV -1 . Here x = sin 8, i?^ = i?e A ", and T p>oc = T p e~ Aa , where T p is the 
temperature at the pole. The general relativistic effects of light deflection are taken into account 
through the ray-tracing function 6 = 6(x), as described in Page (1995). The function lo(0,tp) sets 
the normalization of the flux at each point (6, <fi) of the NS surface 

ft(M) = o 4 Tfl g ii (co^r, (4) 

3 cos z dp + 1 

where, following Heyl & Hernquist (1998b), we have assumed a further dependence of the flux 
on B° . The local temperature on the star is given by T„{d,<j>) = T p [/ (#,</>)] 1/4 , while for the 
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photoabsorption cross section cr(E) we use the analytical fit provided by Morrison & McCammon 
(1988). 

Let now E m \ n and E max define the boundaries of the band of observation, and let us define 
F(Nu) = dEf(E; Nn). Then, if the column density to the source is Ah, the pulse fraction 

in that given band is 

F max (N u ) - F min (AT H ) 

In Figure 1, we plot this function for fixed values of T p = T e ff(ip = 0) and n, while exploring 
the dependence on R/R s . Our results agree with those of DeDeo, Psaltis & Narayan (2000) in the 
limit of Nn — > 0. However, for large column densities, our plots show that the pulsation amplitudes 
can significantly increase. Therefore, neglecting this effect can lead to misleading interpretations. 

In Figure 2, we consider a typical range of column densities inferred for AXPs and show, for 
an intermediate value of R/R s , the dependence of pulse fractions on temperature and degree of 
beaming. We have considered a range of temperatures that is typical for AXPs, if one accounts 
for the fact that processing through the atmosphere typically results in an overestimate of the 
effective temperature by a factor of ~ 2 — 3. Note that the pulsation amplitudes in a given energy 
band increase as T p decreases. This is due to the fact that, as T p decreases, the spectrum is 
shifted towards lower energies. This in turn leads, for the same amount of absorption, to a larger 
suppression of softer photons, and hence to an increase of the pulse fractions. Figure 2 clearly 
shows that a thermal model for highly magnetized NSs can produce pulsation amplitudes which 
are consistent with those observed so far (see Chakrabarty et al. 2000 for a summary). 

Another property that must be accounted for by any viable model for AXPs is the inferred 
size of the emitting region. To estimate the effective emitting areas that an observer would deduce 
for an highly magnetized, cooling neutron star as in our model, we generated a spectrum using 
Equation (|3|). We then used the software XSPEC to convolve it with the instrument response (we 
used CHANDRA) and generate a set of simulated data. These were then fitted using a blackbody 
function. For a star of radius R = {2,3, 4} R s (corresponding to R = {8.25, 12.37, 16.4} km), we 
found, respectively^, Roc ~ {1.3,2.3,3.5} km, significantly smaller than the stellar radius, but 
comparable to the parameters inferred for AXPs. 



3. Conclusions 

We have examined the properties of a model in which the X-ray emission of AXPs is thermal 
radiation from highly magnetized cooling neutron stars. We have shown that the amplitudes of 
the pulsations in a given energy band are an increasing function of the column density to the 
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source. If the typically high column densities inferred for AXPs are taken into account, then even 
a moderate beaming can lead to pulse fractions comparable to those observed so far. 

We have used the software XSPEC to generate simulated data from our model and then used 
these data to find the best-fit blackbody spectrum. We find that the best fit implies an effective 
radius of the star which is typically ~ 5 — 6 times smaller than the actual radius. This is again in 
agreement with the radii of the emitting regions inferred from the observations of AXPs. 

Our results have clear implications for efforts to discern the nature of a variety of X-ray point 
sources. AXPs have variously been interpreted either as accreting sources, or as isolated, strongly 
magnetized neutron stars. If interstellar absorption is not taken into account, then the large 
observed pulse fractions appear to be inconsistent with those predicted by thermal cooling models, 
as found in recent work by DeDeo et al. (2000). However, as we have demonstrated, interstellar 
absorption is able to boost the pulsed fraction measured by a distant observer relative to that 
which would be seen otherwise to the point that a thermal cooling model can easily accommodate 
the pulsed fractions associated with observed AXPs, if the high inferred column densities to these 
sources are properly taken into account. Consequently, we do not believe that pulsed fractions 
discriminate between thermal models for AXP emission from those on which the X-ray luminosity 
is accretion-powered. 

Likewise, interstellar absorption complicates efforts to constrain the sizes of thermally emitting 
areas on the surfaces of neutron stars. The spectrum of the X-ray point source discovered at the 
center of the Cas A supernova remnant (Tannanbaum 1999) has now been discussed separately by 
a number of groups (Pavlov et al. 1999; Umeda et al. 1999; Chakrabarty et al. 2000). Blackbody 
fits to the observed spectra yield emitting areas that are only a small fraction of a neutron star 
surface. These results cannot be explained by other models which interpret the spectra using 
hydrogen atmospheres, but which ignore both the non-uniform temperature distribution expected 
for a strongly magnetized neutron star and the consequences of interstellar absorption. As we 
have indicated, when all these various effects are accounted for, there does not appear to be a 
discrepancy between the effective emitting area of a strongly magnetized neutron star and the 
values which have been derived for objects like Cas A and the AXPs. 



REFERENCES 



Alpar, M. A. 1999, |astro-ph/991222l 



Alpar, M. A. 2000, |astro-ph/0005211 



Chakrabarty, D., Pivovaroff, M.J., Hernquist, L., Heyl, J.S. & Narayan, R. 2000, astro-ph/0001026 
Chatterjee, P., Hernquist, L. &: Narayan, R. 2000, ApJ, in press 
DeDeo, S., Psaltis, D. & Narayan, R., ApJ submitted, astro-ph / 0004266| 
Duncan, R.C. & Thompson, C. 1992, ApJ, 392, L9 



- 6 - 



Ghosh, P., Angelini, L. & White, N.E. 1997, 478, 713 
Greenstein, G., & Hartke, G. J. 1983, ApJ, 271, 283 
Heyl, J. S. & Hernquist, L. 1997a, ApJ, 491, L95 
Heyl, J. S. & Hernquist, L. 1997b, ApJ, 489, L67 
Heyl, J. S. & Hernquist, L. 1998a, MNRAS, 298, L17 
Heyl, J. S. & Hernquist, L. 1998b, MNRAS, 300, 599 
Heyl, J. S. & Kulkarni, S. R. 1998, ApJ, 506, L61 

Hulleman, F., van Kerkwijk, M. H., Verbunt, F. W., & Kulkarni, S. R. 2000, |astro-ph/0002474 

Marsden, D. Lingenfelter, R. E., Rothschild, R. E. & Higdon, J. C. 1999, ApJ, submitted 
( |astFo-ph/9912207D 

Mereghetti, S. & Stella, L. 1995, ApJ, 442, L17 

Morrison, R., & McCammon, D. 1983, ApJ, 270, 119 

Page, D. 1995, ApJ, 442, 273 

Page, D., Sarmiento, A. 1996, ApJ, 473, 1067 

Pavlov, G.G., Zavlin, V.E., Aschenbach, B., Truemper, J. & Sanwal, D. 1999, |astro-ph/9912024 



Perna. R., Hernquist, L., Narayan, R. 2000, ApJ in press, astro-ph/9912297 
Tannanbaum, H. 1999, IAU Circular No. 7246 
Thompson, C. & Duncan, R. C. 1996, ApJ, 473, 322 

Umeda, H., Nomoto, K., Tsuruta, S. & Mineshige, S. 1999, |astro-ph/9910113| 
van Paradijs, J., Taam, R. E. & van den Heuvel, E. P. J. 1995, A&A, 299, 41 
Wang, J.C.L. 1997, ApJ, 486, L119 



This preprint was prepared with the AAS macros v4.0. 



- 7- 




Fig. 1. — Dependence of the pulsation amplitudes on the absorption column density to the source, 
for different values of the NS radius. The flux here is assumed to be observed in the (0.1-10) KeV 
band, the temperature at the magnetic pole is T p = 10 6 K, and the degree of beaming is n = 2. It 
can be seen that a large amount of absorption can significantly increase the observed pulse fractions. 



-8- 




Fig. 2. — Pulsation amplitudes in the range of temperatures and column densities corresponding 
to those typically inferred for AXPs. The flux is taken in the (1-10) KeV band, and the NS radius 
is assumed to be R = 3R S . These values of the pulse fractions are consistent with those observed. 



